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^ ' Abstract 
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We present a calculation of subleading l/m& corrections to spectator effects in the ratios 
of beauty hadron lifetimes in heavy-quark expansion. We find that these effects are sizable 
and should be taken into account in systematic analyses of heavy hadron lifetimes. In 
particular, the inclusion of l/m& corrections brings into agreement the theoretical predictions 
and experimental observations of the ratio of lifetimes of A^-baryon and B d meson. We obtain 
T(B u )/r(B d ) = 1.06 ± 0.02, T{B s )/r(B d ) = 1.00 ± 0.01, T(A b )/r(B d ) = 0.86 ± 0.05. 



1 Introduction 



The hierarchy of lifetimes of heavy hadrons can be understood in the heavy-quark expansion 
(HQE), which makes use of the disparity of scales present in the decays of hadrons containing b- 
quarks. HQE predicts the ratios of lifetimes of beauty mesons, which agree with the experimental 
observations well within experimental and theoretical uncertainties. Most recent analyses give [1, 2] 

r{B u )/r{B d ) \ ex = 1.085 ± 0.017, r(B u )/r(B d ) \ th = 1.07 ± 0.03, (1) 
r(B s )/r(B d )\ ex = 0.951 ± 0.038, r(B s )/r(B d )\ th = 1.00 ± 0.02, (2) 

which show evidence of excellent agreement of theoretical predictions and experimental measure- 
ments. This agreement also provides us with some confidence that quark-hadron duality, which 
states that smeared partonic amplitudes can be replaced by the hadronic ones, is expected to 
hold in inclusive decays of heavy flavors. The heavy-quark expansion also provides a qualitative 
understanding of the hierarchy of lifetimes of charmed hadrons. Since the expansion parameter 
there is much larger, only qualitative agreement is expected. 

The only remaining problem appears in the ratios of lifetimes of heavy mesons and baryons. 
For instance, the low experimental value of the ratio t(A&) /T(B d ) has long been a puzzle for the 
theory. The latest compilation of experimental observations [2] suggests that 

r(A b )/r(B d )\ ex = 0.797 ±0.053, (3) 

which is inconsistent with the first theoretical predictions of this ratio, 0.90 — 0.95 [3, 4, 5] at 
the level of approximately 2a. Only recent next-to-leading order (NLO) calculations of perturba- 
tive QCD [6] and l/m b corrections [7] to spectator effects significantly reduced the discrepancy, 
yielding, as reported in [9], 

T(A b )/r(B d )\ th = 0.88 ±0.05, (4) 

which appears in some agreement with experimental data, mainly due to sizable theoretical and 
experimental uncertainties. The problem of t(Aj,) /r(B d ) ratio could reappear again if future 
measurements at Fermilab and CERN would find the mean value to stay the same with error 
bars shrinking. Upcoming Fermilab measurements of A b lifetime could shed more light on the 
experimental side of this issue. 

Here we expand our calculation of subleading contributions to spectator effects in the \jm b 
expansion [7] to study their impact on the ratios of lifetimes of heavy mesons. We include the full 
charm quark- mass dependence, generalizing our result of Ref. [7]. We also discuss the convergence 
of the l/m b expansion in the analysis of spectator effects by estimating next, l/m^-suppressed 
corrections to the spectator effects in factorization. 

This paper is organized as follows. We set up the relevant formalism and argue for the need 
to compute subleading l/m b corrections to spectator effects in Sect. 2. In Sect. 3 we discuss 
the impact of l/m b corrections to lifetime ratios of heavy mesons and baryons and assess the 
convergence of the l/m& expansion. We present the complete results from \jm\ corrections to 
spectator effects. Finally, we present our conclusions in Sect. 4. Some of the relevant formulas are 
shown in the Appendices. 
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2 Formalism 



The inclusive decay rate of a heavy hadron H b is most conveniently computed by employing the 
optical theorem to relate the decay width to the imaginary part of the forward matrix element of 
the transition operator: 

T(H b ) = -^—(H b \T\H b ), T = lmt fd 4 xT{H cS (x)H cS (0)}. (5) 
ZM Hb J 

Here H e g represents an effective AB = 1 Hamiltonian, 

4G F 



HcS = ~7f Vch £ [Ci{M d ifj) + CMQ U 2 d (a*)J + h.c, 

» d'=d,s,u'=u,c 



(6) 



where the four-quark operators Qi and Q2 are given by 

Qt d> = d' Ll ,u' L CL-fbL, Qt d ' = cl!,u' l 3 L >fb L . (7) 

In the heavy-quark limit, the energy release is large, so the correlator in Eq. (5) is dominated by 
short-distance physics. An Operator Product Expansion (OPE) can be constructed for Eq. (5), 
which results in a prediction of decay widths of Eq. (5) as a series of matrix elements of local 
operators of increasing dimension suppressed by powers of \jm h x : 

m b ) = ^ r Y,m%\H b ) = £^<# b i^ AB= V)i^>- (8) 

2M Hb k k m b 

In other words, the calculation of T(H b ) is equivalent to computing the matching coefficients of the 
effective AB = Lagrangian with subsequent computation of its matrix elements. Indeed, at the 
end, the scale dependence of the Wilson coefficients in Eq. (8) should match the scale dependence 
of the computed matrix elements. 

It is customary to make predictions for the ratios of lifetimes (widths), as many theoretical 
uncertainties cancel out in the ratio. In addition, since the differences of lifetimes should come 
from the differences in the "brown mucks" of heavy hadrons, at the leading order in HQE all 
beauty hadrons with light spectators have the same lifetime. 

The difference between meson and baryon lifetimes first occurs at order 1 / m 2 and is essentially 
due to the different structure of mesons and baryons. In other words, no lifetime difference is 
induced among the members of the meson multiplet. The ratio of heavy meson and baryon 
lifetimes receives a shift which amounts to at most 1 — 2%, not sufficient to explain the observed 
pattern of lifetimes. 

The main effect appears at the 1/m 3 level and comes from the set of dimension-six four- 
quark operators. Their contribution is enhanced due to the phase-space factor 16ir 2 , and induces 
corrections of order 16ii 2 (Aqcd/ m b) 3 = 0(5 — 10%). These operators introduce a difference in 
lifetimes for both heavy mesons and baryons. The effects of matrix elements of these operators, 



1 We choose a heavy-quark pole mass as an expansion parameter, which is consistent with most of the previous 
calculations. Of course, final predictions for the lifetime ratios are independent of that choice [8]. 
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commonly referred to as Weak Annihilation (WA), Weak Scattering (WS), and Pauli Interference 
(PI) have been computed in [4, 11, 12, 13, 14, 15, 16] at leading order in perturbative QCD, and, 
more recently, including NLO perturbative QCD corrections in [6] and 1/mj, corrections [7]. 

The contributions of these operators to the lifetime ratios are governed by the matrix elements 
of AB = four-fermion operators 



rj- ryu _j_ rr'd _j_ <~r~S 

-'spec — -'spec ' -'spec ' "'spec) 



(9) 



where the % are 



spec 



-^spec 



G 2 F m 2 \V bc \ 2 (l - z) 



2tt 

G\m\\V hc \\\-zf 



^{ (c\ + 4) 0\ + 2 Cl c 2 0r + 5 u 1/m + 6 u l/m2 }, 



4% 



[l + z)Of + -{l + 2z)Oi 



spec 



+ (N c c\ + 2cic 2 
G 2 F m 2 b \V bc \ 2 VT=Iz~ 



l + z)Of + |(l + 2z)Ctf 



+ °l/m + °l/m? 



(10) 



47T 



or 



2z)O s 2 



+ (N c c\ + 2 Cl c 2 



01 +-(i + 2z)O s 2 



/m 2 



Here the terms 5y m and 5\, m 2 refer to l/rrib and 1/mj; corrections to spectator effects, which we 
discuss below. Note that we include the full z = m 2 /m 2 dependence, which is fully consistent only 
after the inclusion of higher l/m^ corrections. The operators Oi and Oi in Eq. (10) are defined as 



Of = ba»(l - 75)^^(1 - 7s)?j, <5f = kYlbbjqa^l - 75)^. 



(11) 



In addition, there is a dimension-6 operator O3 = 6j7 M (l — 75)^^7^75^ whose contribution is 
always suppressed by powers of the light-quark mass. It gives a negligible contribution to the 
ratios of lifetimes of heavy hadrons (see also Sect. 3). 

In order to assess the impact of these and other operators, parametrizations of their matrix 
elements must be introduced. The meson matrix elements are 



{B q \0\\B q ) = f 2 m 2 Bq (2 Cl + ^) , (B q \0\\B q ) = f 2 Bq m 2 B B u 



(B q \Ol\B q ) = -fl 



m. 



{m h + m q ) V N c 



B-, 



■2, , + — ) 1 - I 2ei + ^ 



B 1 



(12) 



(B q \Ol\B q ) = -flm 



m B 1 

B ~ 2^2 H -^1 

" |_(to 6 + m g ) " 2 



Here the parameters Bi and are usually referred to as "singlet" and "octet" parameters (see Ap- 
pendix A for the explanation). Expressed in terms of these parameters, the lifetime ratios of heavy 
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mesons can be written as 



r(B u )/r(B d ) = 1 + 16tt 2 [G? (m^m,) + G°° (m b )ei(m b ) 

+ G s 2 s (m b )B 2 (m b ) + G?(m b )e 2 (m b )] + 5 1/m , (13) 

where the coefficients G were computed at NLO in [15, 6]. 8\/ m represents spectator corrections 
of order 1 / m b and higher, which we present below in this paper. 

Estimates of the matrix elements of four-fermion operators in baryon decays are not easy. 
Following [4], the parameter B is used to account for the deviation of the A b wave function from 
being totally color-asymmetric, 

(A b \Ol\A b ) = -B(A b \Ol\A b ) = jf Bq m Bq m Ab r, 

(A b \O q 2 \A b ) = -B{A b \6l\A b ) = jf Bq m Bq m Ab 5. (14) 

Note that 5 = 0(l/m b ), which follows from the heavy-quark spin symmetry. It needs to be in- 
cluded as we consider higher-order corrections in l/m b . B — 1 in the valence approximation. While 
these parameters have not been computed model-independently, various quark-model arguments 
suggest that the meson and baryon matrix elements are quite different. Thus a meson-baryon 
lifetime difference can be produced. In general, one can parametrize the meson-baryon lifetime 
ratio as 

r(A b )/r(B d ) ~ 0.98 - {d 1 + d 2 B)r - (d 3 ei + d 4 e 2 ) - {d b B^ + d 6 B 2 ) 

~ 0.98 - m\{d\ + d' 2 B)r - m 2 b [(d^ + d' 4 e 2 ) + (d'^ + d' 6 B 2 )} , (15) 

where in the last line we scaled out the coefficient ml emphasizing the fact that these corrections 
are suppressed by l/mf compared to the leading ml effect. The scale-dependent parameters 
(di(m b ) = {0.023, 0.028, 0/16, -0.16, 0.08, -0.08} at NLO) are defined in [4]. The parameter 

r = ip b g(0) I ip b q(0) is the ratio of the wave functions at the origin of the A b and B q mesons. 

It is interesting to note that in the absence of l/m b corrections to spectator effects, it would 
be equally correct to substitute the 6-quark mass in Eq. (15) with the corresponding meson and 
baryon masses, so 

r(A b )/r(B d ) ~ 0.98 - m 2 ^ + d' 2 B)r - m% d [(dfa + d' A e 2 ) + (d^B, + d' G B 2 )\ , (16) 

which reflects the fact that WS and PI effects occur for the heavy and light quarks initially bound 
in the B d meson and A b baryon, respectively. While correct up to the order 1/mf, these simple 
substitutions reduce the ratio of lifetimes by approximately 3-4%! We take this as an indication 
of the importance of bound-state effects on the spectator corrections, represented by subleading 
l/m b corrections to spectator operators. 
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Figure 1: Kinetic corrections to spectator effects. The operators of Eqs. (17) 
and (19) are obtained by expanding the diagrams in powers of spectator's 
momentum. 

3 Subleading corrections to spectator effects 

We computed the higher order corrections, including charm quark-mass effects, to Eq. (10) in the 
heavy-quark expansion, denoted below as 8\j m and $\i m i- 

3. 1 1/mh corrections 

The 1/mb corrections to the spectator effects are computed, as in Ref. [7], by expanding the 
forward scattering amplitude of Eq. (5) in the light-quark momentum and matching the result 
onto the operators containing derivative insertions (see Fig. 1). The Sy m contributions can be 
written in the following form: 



°1, 



/m 



v ' 1 — z 1 — z 



1 + z 



I III 



si 



I in 



l-z 3 1-z 1 3 V ; V 2 3 



+ [N c c 2 2 + 2cic 2 
16z 2 



52' 



R d o + 



2 1 + z + 10z 2 -* 2 



l-Az RS ° + 3 



+ (N c c 2 2 + 2cic 2 



l-z u 3 
2 1 - 2z + 16z 2 



1 - z 



[IT) 



l-4z 



itf + -(l + 2z) (itf-itf 



1 -4z 



l-4z 
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where the following operators contribute 




1 




1 



2^7^(1 - 75)£> Q Mj7ju(l - l5)D a qj, 



ml 



ml 




1 



m, 



^(1 - 75)6i9j(l - 7s) 



(18) 



Here i?f denote the color-rearranged operators that follow from the expressions for R\ by inter- 
changing color indexes of bi and qj Dirac spinors. Note that the above result contains full QCD 
6-fields, thus there is no immediate power counting available for these operators. The power 
counting becomes manifest at the level of the matrix elements. 

In the z — > limit the result of Eq. (17) reduces to the set of l/m& corrections obtained 
by us in Ref. [7]. We complete the calculation presented there by computing the contribution 
of Eq. (17) to the ratios of heavy- meson lifetimes. As before, we use factorization to guide our 
parametrizations of A b and meson matrix elements, but, contrary to Ref. [7], we elect to keep 
corrections that would vanish in the vacuum-insertion approximation. Our parametrizations for 
meson and baryon matrix elements are presented in Appendix B. Inserting Eqs. (B.2) and (B.5) 
of Appendix B into Eq. (17) we obtain an estimate of the effect. We postpone the presentation of 
the answer, complete with our estimate of its uncertainties, until the end of this section, after we 
assess the convergence of the l/m b expansion for the spectator corrections. 

Numerically, the set of l/m b corrections does not markedly affect the ratios of meson lifetimes, 
changing the t(B u )/t(B ( i) and the t{B s )/t{B ( ]) ratios by less than half a percent. The effect 
is more pronounced in the ratio of A b and B^ lifetimes, where it constitutes a 40 — 45% of the 
leading spectator contribution represented by WS and PI effects (see Eqs. (15,16)), or an overall 
correction of about —(3 — 4)% to the r(A b )/r(B d ) ratio. While such a sizable effect is surprising, 
the main source of such a large correction can be readily identified. While the individual l/m b 
corrections to the matrix elements representing WS and PI are of order 20%, as expected from the 
naive power counting, they contribute to the A b lifetime with the same (negative) sign, instead 
of destructively interfering as do the contributions to Eqs. (15,16) representing WS and PI [5], as 
explained in [7]. This conspiracy of two ~ 20% effects produces a sizable shift in the ratio of the 
A b and 5-meson lifetimes. 

3. 2 1/^5 corrections 

In order to see how well the l/m b expansion converges in the calculation of lifetime ratios, we 
compute a set of Si/ m 2 corrections to spectator effects. It is expected that at this order even 
more operators will contribute, so the exact prediction of lifetime differences becomes unfeasible. 
Therefore, we shall use our calculation only to see to what extent we can trust our results for 
the l/m b corrections. We will parametrize \jm\ corrections similar to our parametrization of 
l/m b effects above and will use the factorization approximation to assess their contributions to 
the lifetime ratios. 

There are two classes of corrections that arise at this order. One class involves kinetic correc- 
tions which can be computed in a way analogous to the previous case by expanding the forward 
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scattering amplitudes in the powers of spectator momentum (see Fig. 1). A second class involves 
corrections arising from an interaction with background gluon fields. The kinetic corrections can 
be written as: 



mil — z 



l-z) 



+ 2cic 2 



mi 1 — z 



l-z) 



mll + z + 2z 2 



mi 



l-z 



ml Az 2 ml 2(1 + 2z) 
Of H %■- 0% -\ '-Oi 



m 2 ! — z 



mi 



Az 2 (7z-5), , 8z 2 (Az-3) 



(i- z y 



■W? + 



+ (N c c 2 2 + Cl c 2 



(l-z) 2 
ml 1 + z + 2z 2 



W. 



f + ^- (wf - wf) 

l — z y ' 



of + 



mi 



rn 



% ^ , rn 2 , 2(1 + 2z) 



m 2 1 — z 



O a , + 



mi 



Of 



Az 2 (7z-b)~ d , 8z 2 (Az-?>)~ d , 8z 2 ( ~ d , 



°l/m 2 



= Ci 



(1-*)= 

TO^ 1 — 2z 
ml I — Az 



(l-zf 



, m 2 sl 8z 2 , m*, 2(1 + 2z), 

1 + -Ti + 2 o 3 



m% 1- Az 



mi 



8z 2 (lQz-5) 
+ (1-Az) 2 

+ (N c c 2 2 + 2cic 2 

8z 2 (lQz-h) 



Wl + 



lQz 2 (10z-3) 
(1-Az) 2 



W° + 



3 

IQz 2 
1-Az 



(wf - Wi) 



m 2 s , l — 2z 
mil — Az 



0\ + 



m s , oz 
ml 1 — Az 



O* + 



ml 2(1 + 2z) 



s , 



mi 



(l-Azf 



16z 2 (l(k-3)~ , lQz 2 ,~. ~,x 
W i H 7T a \i W 2 + i — (W3 - W4 J 



1-4* 



(19) 



We again retain the dependence on quark masses in the above expression, including the terms 
proportional to light-quark masses (and thus multiplied by the operators Oj). The operators in 
Eq. (19) are defined as 

W? = \hY(l ~ ^)D a D\q jltl (l - lb )D a D pqj , W% = \b a ^D a D%q n ,(1 - lz )D a D mj , 



ml 



ml 



W% = — & i7 "(l - -y 5 )D"D a b i q jlv (l - K)D»D aqj , W\ = -^(1 - <y 5 )D a b t q 3 (l - 75) (20) 
m b mj, 

where, as before, Wf denote the color-rearranged operators that follow from the expressions for 
Wf by interchanging the color indexes of foj and qj Dirac spinors. The parametrization of matrix 
elements of these operators can be found in Appendix B. 

In addition to the set of kinetic corrections considered above, the effects of the interactions of 
the intermediate quarks with background gluon fields should also be included at this order. The 
contribution of those operators can be computed from the diagrams of Fig. 2, resulting in 



7- 

-^spcc,G 



rru 1 rrd' I rrs' 

spec,G ' spec,G ' spcc,G ' 



(21) 
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where 




Figure 2: Corrections to spectator contribution coming from the interactions 
with background gluon fields. 



spec,G 



rr-d! 

spec,G 



spec,G 



47T 



[c\ [(1 - z 2 )Pf - (1 - z 2 )P( + 2z{\ - z)Pf + Az 2 Pf 



+2 Cl c 2 z [(1 - z)Pf + (1 - z)Pf + 2zPf + 2zPf] } , 



+4cic 2 z 



PS + Pi + P° 



P S 



} 



The local four-quark operators in the above formulas are shown in Eq. (25). 



P? 



P q 



pi 
^3 



pi 



pi 



pi 



m 



^(1 + 75)^^^-^(1-75)%, 
kY(l+l 5 )h ^,^(1 -75)%, 



Pr q = 


m b 


-75)?/, 


n = 


-^b t Y(l ~ ls)D v D«b k t a kl G; u q lla (l 
m b 


- 7s)?i, 


p$ = 






P?o = 


h^h t^G^Yii - 75 ) ft . 





Analogously to the previous section, we parametrize the matrix elements as 

(B q \Pf\B q ) = \f B m\a\ " l) , (A b |^|A 6 > = \f Bq m Bq m Kb at - l) . (26) 

We set a* = a A = 1 GeV 2 to obtain a numerical estimate of this effect. It is clear that no 
precise prediction is possible with so many operators contributing to the ratios of lifetimes. This, 
of course, is expected, as the number of contributing operators always increase significantly with 
each order in OPE. Fortunately, at least in factorization, the effect of 1 /m 2 corrections is tiny, a 
fraction of a percent. The contribution of l/m 2 -suppressed effect is included into the coefficient 
ko defined in the next section. We can evaluate the contribution of both l/mt, and 1/m 2 by 
performing two analyses. First, we can study the perturbative scale dependence of the final result 
using "factorized" values for the matrix elements 2 . We shall see that the scale dependence of the 
resulting effect is very mild. Second, we shall randomly vary the parameters describing the matrix 
elements by ±30% around their "factorized" values to see the uncertainty of our result. 



3. 3 Phenomenology of heavy hadron lifetimes 

Let us now discuss the phenomenological implications of the results presented in the previous 
sections. As usual in the OPE-based calculation, next-order corrections bring new unknown 
coefficients. In our numerical results we assume the value of the 6-quark pole mass to be = 
4.8 ± 0.1 GeV and f B = 200 ± 25 MeV. 

The resulting lifetime ratios are parametrized according to Appendix B. For instance, in the 
case of T(A b )/r(B d ) we can write 

= k + ]T k iPi , (27) 

h.o. i=l 

where pi are the parameters defined in Appendix B, 

p t = {i? 1 ,E 2 ,6 1 ,e 2 ,r,^,Er,^,/3 ,^ o, fa, fa, fa, fa, fa, fa, ^, ft A , ^ A , & A , & A , /5 3 A , ^ , (28) 

2 In the absence of model-independent calculations of the matrix elements described above, therejs a number of 
ways to estimate them. We chose to use the quark model to relate the matrix elements of Ri and Ri operators to 
the matrix elements of dimension 6 operators Oi, which are parametrized in terms of Bi and ej. We then use lattice 
QCD predictions for Bi and Ci to obtain numerical estimates of our corrections. We shall refer to this procedure 
as "factorization." 



r(B d ) 
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Coefficients fc, for the parameters listed below. All entries 


are x 10 


-3 






Bi 


B 2 


ei 


£2 


r 


(5 


Br 


BS 


0o 


A> 


0i 


0i 


m b /2 


-7.82 


9.19 


-158 


183 


-27 '.4 


3.76 


-33.3 


11.1 


0.001 


1.95 


-0.014 


-24.0 


nib 


-7.96 


9.31 


-160 


196 


-22.4 


3.39 


-27.8 


12.4 


0.031 


1.78 


-0.382 


-21.9 


2m b 


-7.45 


8.63 


-160 


201 


-18.9 


3.33 


-24.6 


13.0 


0.079 


1.69 


-0.976 


-20.9 




02 


02 


03 


03 




/?o A 














m b /2 


-0.002 


-3.68 


0.002 


3.68 


-0.144 


-0.048 


-12.7 


-6.85 


-0.818 


-0.272 


-0.818 


-0.272 


m b 


-0.058 


-3.36 


0.058 


3.36 


-0.132 


-0.039 


-11.2 


-4.59 


-0.747 


-0.223 


-0.747 


-0.223 


2uib 


-0.150 


-3.20 


0.150 


3.20 


-0.125 


-0.030 


-10.4 


-3.09 


-0.711 


-0.170 


-0.711 


-0.170 



Table 1: Coefficients fcj appearing in Eq. (27) for r(A b )/r(S d ). 





Coefficients fcj for the parameters listed below. All entries arc xlO 3 . 
B x B 2 £l e 2 0o 0o 0i A 02 02 03 03 


rrib/2 
m b 
2nib 


82.1 18.8 -782 206 0.001 1.95 -17.6 124 -0.002 -3.68 0.002 3.68 
60.8 13.3 -727 211 0.031 1.78 -6.58 107 -0.058 -3.36 0.058 3.36 
46.5 9.98 -693 213 0.079 1.69 0.661 98.9 -0.150 -3.20 0.150 3.20 



Table 2: Same as Table 1 for T(B u )/r(B d ). 

and hi are their coefficients, ko = 1.00 in the case of mesons, while for A b ko = {0.991, 0.981, 0.975} 
for (j, = {m b /2, m b , 2 m b }, respectively. They appear in Table 1, and a similar procedure is followed 
for T{B u )/r{Bd) and t{B s )/t{B (J ), whose coefficients are shown in Tables 2 and 3, respectively. 

Since all three heavy mesons belong to the same SU(3) triplet, their lifetimes are the same 
at order l/mf. The computation of the ratios of heavy meson lifetimes is equivalent to the 
computation of £/-spin or isospin- violating corrections. Both 1 /m^-suppressed spectator effects 
and our corrections computed in the previous sections arise from the spectator interactions and 
thus provide a source of [/-spin or isospin-symmetry breaking. We shall, however, assume that 
the matrix elements of both l/mf and \jm\ operators respect isospin. The ratio of lifetimes of 
B s and B d mesons involves a breaking of [/-spin symmetry, so the matrix elements of dimension-6 
operators could differ by about 30%. We partially take this breaking into account by retaining 
light quark-mass dependence. Nevertheless, we still introduce different B— and e-parameters to 
describe B s and B d lifetimes. We will also assume that, apart from the explicit light quark-mass 
effects, matrix elements of 1/mf operators respect [/-spin symmetry. 

Since the computation of 1/ m b corrections to spectator effects was performed at leading order 
in QCD, we also studied the scale dependence of the result. A comparison between the results 
computed with leading logarithmic accuracy (LL), next-to- leading log accuracy (NLL) and the 
total NLL including l/m b and I /ml corrections is plotted for r(A b )/r(B d ) in Fig. 3 for a renor- 
malization scale parameter interval m b /2 < /i < 2m b , assuming "factorization", i.e. setting all the 
/3-parameters in Eqs. (B.2), (B.5) to the value calculated using our "factorization" ansatz. Figs. 4 
and 5 show analogous plots for r{B u )/r{B ( j) and t(B s ) /r{B d ). As can be seen from these graphs, 
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Coefficients fc, for the parameters listed below. All entries 


are x 10 3 . 




B l 


B 2 


£i 


£2 


B{ 


B s 2 




e 2 


m b /2 


-7.82 


9.19 


-158 


184 


11.2 


-12.4 


180 


-230 


mb 


-7.96 


9.31 


-160 


196 


13.5 


-16.8 


182 


-243 


2rrib 


-7.45 


8.63 


-160 


201 


14.7 


-19.1 


183 


-249 




Po 


Po 


Pi 


Pi 


Pi 


P2 


Ps 


Ps 


m b /2 


-0.003 


-5.72 


0.011 


18.9 


0.0013 


2.20 


-0.0013 


-2.20 


m b 


-0.091 


-5.23 


0.300 


17.3 


0.0349 


2.00 


-0.0349 


-2.00 


2mb 


-0.233 


-4.97 


0.768 


16.4 


0.0892 


1.91 


-0.0892 


-1.91 



Table 3: Same as Table 1 for T(B s )/r(B d ). 
the scale dependence of the predictions is very mild. Next, we fix the scale in our calculations 




4 6 



Ijl (GeV) 

Figure 3: LL (dashed line), NLL (dash-dotted line), and NLL+l/mf, (solid 
line) contributions to t(A&) /r{Bd) are plotted vs. the scale parameter \i in 
GeV. 

and vary the values of parameters of the matrix elements. We adopt the statistical approach for 
presenting our results and generate 20000-point probability distributions of the ratios of lifetimes 
obtained by randomly varying our parameters within a ±30% interval around their "factoriza- 
tion" values, for three different scales [i. The decay constants Jb and b-quark pole mass are 
taken to vary within la interval indicated above. The results are presented in Figs. 6, 7, and 8. 



12 



1.100 
1.075 
1.050 



L 

\ 

\ 
\ 

\ 



1.025 - 



1.000 - 



0.975 - 



-I I I I I I I I I I i_ 



-I I 1- 



6 



8 



/x (GeV) 



Figure 4: Same as Fig. 3 for T(B u )/r(B d ). 
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Figure 5: Same as Fig. 3 for t{B s )/t{B^). 
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These figures represent the main result of this paper. They show a nice agreement between the 
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0.75 0.80 0.85 0.90 0.95 0.75 0.80 0.85 0.90 0.95 0.75 0.80 0.85 0.90 0.95 

r(A b )/r(B d ) 

Figure 6: Histograms showing the random distributions around the central 
values of the /n , m b , B, 5, e, Pi and /3« parameters of Eqs. (B.2) and (B.5) 
contributing to r(Aft)/r(5 £ j). Three histograms are shown for the scales /J = 
mb/2 (a), fi = nit, (b), and /x = 2mj (c). The shaded curves, presented for 
the convenience, represent current experimental result. Their normalization is 
arbitrary. 




1.04 1.06 1.08 1.10 1.12 1.02 1.04 1.06 1.08 1.10 1.02 1.04 1.06 1.08 1.10 

T(B u )/r(B d ) 

Figure 7: Same as Fig. 6 for r(B u ) /r{Bd). 

experimental results and theoretical predictions for all of the ratios of lifetimes of heavy hadrons. 
In fact, it is a relatively "large" ratio of t(B s )/t(B c i) that can pose a potential problem if future 
measurements at Fermilab and CERN would find the mean value to stay the same with its error 
bars shrinking. Upcoming Fermilab and CERN measurements would again clarify this issue. 

It also appears that the theoretical errors in the ratio of t(B s )/t(Bci) are much smaller than 
the other two. We believe that the main reason for it is the fact that l/m\ and higher specta- 
tor corrections to both B s and B& inclusive widths are given by the same type of diagrams. In 
other words, apart from the different charm-quark mass dependence (which comes at the level of 
m 2 c jm\ corrections to 1/mf effects), both widths of B s and Bd are affected by the weak scatter- 
ing contributions in the same way. Thus, besides possible S"?7(3)-breaking effects in the matrix 
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0.94 0.96 0.98 1.00 1.02 1.04 0.94 0.96 0.98 1.00 1.02 1.04 0.94 0.96 0.98 1.00 1.02 1.04 



r(B s )/r(B d ) 

Figure 8: Same as Fig. 6 for r(B s )/r(B d ). 

elements of the 4-fermion operators, the contributions of spectator interactions to the numerator 
and denominator of the ratio of B s and widths cancel each other to a large degree. This effect 
persists for the class of l/m& corrections to spectator effects discussed here as well. 

4 Conclusions 

We computed subleading l/m& and 1 / ml corrections to the spectator effects driving the difference 
in the lifetimes of heavy mesons and baryons. Thanks to the same 167T 2 phase-space enhancement 
as 1/m^-suppressed spectator effects, these corrections constitute the most important set of 1/mf- 
suppressed corrections. 

We showed that a set of l/m^-corrections to spectators effects can be parametrized by eight 
new parameters $ and (i = 1, 4) for heavy mesons and eight new parameters /3 4 A and for 
heavy baryons. Although model-independent values of these parameters will not be known until 
dedicated lattice simulations are performed, we present an estimate of these parameters based on 
quark model arguments. Compared to [7] , we kept all matrix elements, including those vanishing 
in the factorization approximation. This is important, as the Wilson coefficients of those operators 
are larger than the ones multiplying the operators whose matrix elements survive in the N c —>■ oo 
limit. We also perform studies of convergence of this expansion by computing a set of 
corrections to spectator effects and estimating their size in factorization. The expansion appears 
to be well- convergent. 

The main result of this paper are Figs. 6, 7, and 8, which represent the effects of subleading 
spectator effects on the ratios of lifetimes of heavy mesons and baryons. We see that subleading 
corrections to spectator effects affect the ratio of heavy meson lifetimes only modestly, at the 
level of a fraction of a percent. On the other hand, the effect on the Ab-Bd lifetime ratio is quite 
substantial, at the level of —3%. This can be explained by the partial cancellation of WS and PI 
effects in baryon and constructive interference of 1/mj, corrections to the spectator effects. 

There is no theoretically-consistent way to translate the histograms of Figs. 6, 7, and 8 into 
numerical predictions for the lifetime ratios. As a useful estimate it is possible to fit the histograms 
to Gaussian distributions and extract theoretical predictions for the mean values and deviations 
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of the ratios of lifetimes. Predictions obtained this way should be treated with it is not 

expected that the theoretical predictions are distributed according to the Gaussian distribution. 
This being said, we proceed by fitting the distributions to Gaussian and, correcting for the scale 
uncertainty [varied in the interval (rrib/2, 2to&)] by inflating the error bars, extract the ratios 
T(B u )/r(B d ) = 1.06 ± 0.02, T(B s )/T(B d ) = 1.00 ± 0.01, and T(A b )/r(B d ) = 0.86 ± 0.05. In these 
equations the mean values are taken at the scale /i — m^. This brings the experimental and 
theoretical ratios of baryon and meson lifetimes into agreement. 
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Appendix A. Four-quark operators 

This paper uses the basis of 1/mf operators consistent with Ref. [3], but different from Refs. [4, 6], 
given the relative simplicity of the higher-order operators in our basis. In the above calculations 
we used the following notation for local four-quark operators: 

0\ = kY(l ~ 75)^7^(1 - 7s)?j, 0\ = biYlbhqjl^ ~ 7s)?j, 

O q = 1^(1 - -f 5 )biq jWqj , (A.l) 

as well as Of operators that follow from the expressions for Of by interchanging the color indexes 
of hi and qj Dirac spinors. Refs. [4, 6] use more familiar "color-singlet" and "color-octet" operators 

Q q = b L -f^q L q L ^b L , Q q s = b R q L q L b R , 

T q = b Ll ,t a q L q Ll H a b L , Tj = b R t a q L q L t a b R . (A.2) 
These operators can be expressed in our basis as 

Q q = \oi Ql = -\o\ - \oi 

T q = -^-d\ + \o\, T$ = (\d\ + O q ) - I (\o\ + Ol) . (A.3) 
8N C 1 8 1 b 8N C V2 1 V 8 V2 1 V y J 

Note that in the predictions for lifetime ratios, large matrix elements of the "color-singlet" op- 
erators Q q and Q q s are multiplied by a relatively small combination of Wilson coefficients, while 
matrix elements of "color-octet" operators come with large Wilson coefficients [4, 14]. Neverthe- 
less, the interpretation of matrix elements obtained in the vacuum insertion approximation (i.e. 
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representing the matrix elements of four-fermion operators as products of two matrix elements of 
current operators separated by a vacuum state) is relatively clear in this basis. For the dimension-6 
operators this procedure introduces four new scale-dependent parameters -Bj(/x) and €i(fj), 

= ^.GO, j^WW = f ^M- (a.4) 

These parameters can be computed using QCD sum rules, quark models, or on the lattice. A 
compilation of various estimates of these parameters can be found in [17]. 



Appendix B. Estimate of matrix elements 



We parametrize relevant matrix elements using factorization results as a guiding principle (a 
similar approach was used in [4, 7, 18]). Of course, these matrix elements would eventually have 
to be computed model-independently, using lattice QCD or QCD sum rules. In the absence of 
model-independent results, we obtain numerical predictions for the ratios of lifetimes, by relating 
matrix elements of the computed operators constituting l/m& and 1/mjj effects to the matrix 
elements of the operators of dimension 6 (Oi )2 ) using quark model. We then use lattice QCD 
predictions for the parameters Bi and q and look at the ratio of the subleading- and leading-order 
spectator effects. We believe that some errors introduced by our use of the quark model should 
partially cancel in these ratios. The leading-order meson matrix elements are 



(B q \0\\B q ) 
(B q \O q 2 \B q ) 

(B q \Ol\B q ) 



fiy Bg ( 2 Cl + ^) , (B q \dl\B q ) = flmlB,, 



r-2 2 

-fB a ™ Ba 



ffl 



rl 2 
~fB„ m B, 



(m b + m q y 



2e 2 + 



B? 



+ « 2d + 



Bi 



(B.l) 



-^B 2 + -B 1 
(m b + m q ) 



1 

2" 



(B q \O q 3 \B q ) = (B q \O q 2 \B q ), (B q \O q 3 \B q ) = (B q \O q 2 \B q ). 

The parameters Bi and represent matrix elements of "color-singlet" and "color-octet" operators. 
The above parameters can be computed in QCD sum rules, quark models, or on the lattice. 
Naively, one expects that in the large- N c limit Bi ~ 0(1), t{ ~ 0(1/N C ). 

A parametrization of l/m 6 -suppressed corrections could, in principle, be quite simple: after 
scaling out the proper m^-dependence each new matrix element brings in a new parameter. This 
parametrization could be useful for lattice studies of these matrix elements. Here we shall choose 
a slightly different approach. Using the quark model to guide our parametrization, we relate our 
matrix elements to the matrix elements of singlet and octet operators using Fierz relations. We 
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find 



K\B q ) 
Rl\B q ) 
R\\B q ) 
R\\B q ) 

m\B q ) 

R q 2 \B q ) 



Rl\B q ) 



Ri\B q ) 




^fl q m l q Po 



m, 



mi 



- 1 



2 

\fl q ™ 2 B q Pl 



m 



Bo 



mi 



1 

12' 




(B.2) 



- 1 



1 f 2 2 n ( m l q / 



1/2 m | (& +2/ 3. 



B„ 




- 1 



In the large N c - limit ~ and ~ 

Pa = + Si/2), A = /3 2)3 = /Kg 



0(1/N C ). We can use the quark model to relate 

The same 



- 4B 2 )(1 + m 2 /m| 9 ) + + 2B 2 ) 

relations hold for $ with i?j — > e^. Here (3[ are used to parameterize the ratios of the matrix 
elements of the four fermion operators and 1 / m^-suppressed operators. In our numerical studies 
we took the parameters B{, q from the lattice calculations of matrix elements of dimension-6 
operators [6] and varied f3[ around its quark-model value of unity by ±30%. 

A similar logic can be used to estimate l/m 2 -suppressed matrix elements. We simplify the 
result by using the factorization approximation, i.e. neglecting all the matrix elements of octet 
operators, which should suffice for our purposes of studying the convergence of the series. We 
obtain 



{B q m\B q ) 
{B q m\B q ) 

(B q \W^\B q ) 



fl q m l q ( 




4iV c ^ 


m 2 b 




m k 




m 2 b 






8N C 
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8 


V m l 


fk< ( 




8N C { 


m 2 b 




m % 




m 2 b 



(B q \W«\B q ) = -f - 1 , (B.3) 



(B q \W!\B q ) = 



(B q m\B q ) 



(B q \W!\B q ) = (B q \W§\B q ), (B q \W!\B q ) = (B q \W§\B q ). 

In a similar way, we used the quark-diquark valence-quark model to guide our parametrizations 
of baryon matrix elements. The leading order matrix elements are [4, 7] 



(A b \Ol\A b ) 


= -B(A b \Ol\A b ) 


B f2 




(A b \O q 2 \A b ) 


= -B{A b \O q 2 \A b ) 


= ^fB q m B q m Kb 5, 


(B.4) 


(MOl\A b ) 


= -B(A b \d q 3 \A b ) 


= (A*|Ol|A 6 >. 





where r 



■06^(0) / 'tPbqW * s the ratio of the wave functions at the origin of the A b and 

B q mesons, and B — 1 in the valence-quark model. Estimates of r vary from 0.1 to 1.8 and 
can potentially be larger [4]. Note that 5 = 0(l/m b ), which follows from the heavy-quark spin 
symmetry. We use the same principles to parametrize the l/m fe -suppressed contributions, 



(M\R q o\A b ) 


= -f^(A b \Rl\A b ) = 
Po 


00 i-2 

-^f Bq m B q m Ab 


( m l b 
\ ™ 2 b 


(A b \Rf\A b ) 


= -S.(A b \Rl\A b ) = 
Pi 


-^fB q mB q m Ab 


( m l b 

V m l 


(A b \R q 2 \A b ) 


= -^(A b \Rl\A b ) = 

P2 


P2 r2 

-^fB q mB q m Ab 


( m l b 

\ m b 


(Ab\Rt\A b ) 


= -^(A b \Rl\A b ) = 

P3 


(3 A ( 
-^fl q m Bq m Ab ^ 


m l b 

m 2 b 



- 1 



(B.5) 



Note that the matrix elements contain the parameters (3$ and /3q, which are of the order of l/m b . 
Thus, these matrix elements only contribute at the order \jm\ corrections to spectator effects (see 
below). In the valence-quark approximation these parameters are proportional to 5. Analogously 
to the meson case, the quark model relates $ = (3' A B5, = $$5, p\ A = (3[ A Br, (3 A = (3 A 'r, 
P*iz = & [(r + 4 5/3)(l + ml/m\ b ) - 2 5/3] , & A 3 = [(r + 4 J/3)(l + mg/m|J - 2 5/3] . We 
use the valence-quark approximation for l/mf corrections. We obtain 



(A b \W?\A b ) = -<A 6 |W^|A 6 ) 



fl q m Bq m Ab 



96 



ml, 



mi 
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(A b \W«\A b ) = -{A b \W$\A b ) = (B.6) 

(A b \W!\A b ) = -{A b \W!\A b ) = -{A b \W$\A b ), 
where we set all the terms of order 0(l/m^) in the matrix elements to zero. 
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